Magnetic order in the pseudogap phase of high- Tq superconductors 
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One of the leading issues in high-Tc superconductors is the origin of the pseudogap phase in 
underdoped cuprates. Using polarized elastic neutron diffraction, we identify a novel magnetic 
order in the YBa2Cu306+a: system. The observed magnetic order preserves translational symmetry 
as proposed for orbital moments in the circulating current theory of the pseudogap state. To date, 
it is the first direct evidence of an hidden order parameter characterizing the pseudogap phase in 
high-Tc cuprates. 



In optimally and underdoped regimes, high-Tc C0 P_~ 
per oxides superconductors exhibit a pseudogap statellj, 
LH with anomalous magnetic^], transport p, 
thermodynamic^] and optical^] properties below a tem- 
perature, T*, large compared to the superconducting 
transition temperature, T c . The origin of the pseudo- 
gap is a challenging issue as it might eventually lead 
to identify the superconducting mechanism [[J. Two ma- 
jor classes of theoretical models attempt to describe the 
pseudogap state: in a first case, it represents a precursor 
of the superconducting ei-wave gap@, with preformed 
pairs below T* which would acquire phase coherence be- 
low T c 0,^3- I n a second approach, the pseud ogap is 
associated either with an ordered 0, 0, 0, Q, Hilllrj] 
or a disordered phase Q, 0, 0] competing with the 
SC one. The order parameter, associated with these 
competing phases may involve charge and spin den- 
sity waves [l^ llnL llfij or charge currents flowing around 
the CuO? square lattice, such as D-charge density wave 
(DDW) O or orbital circulating currents (CC) [HE2]. 

Most of these phases break the translation symmetry 
of the lattice (TSL). Therefore, they may induce charge, 
nuclear or magnetic superstructures that can be probed 
by neutron or X-ray diffraction techniques. In contrast, 
CC phases [H E2 preserve the TSL as they correspond 
to 4 or 2 current loops per unit cell (referred as 0/ and 
0/j phases, respectively). These charge currents could 
be identified by virtue of the pattern of ordered orbital 
magnetic moments pointing perpendicularly to the Cu02 
planes. These orbital magnetic moments should be de- 
tectable by neutron diffraction. Although the TSL is 
preserved, the magnetic signature of the CC phase does 
not reduce to ferromagnetism: the loops are staggered 
within each unit cell corresponding to a zero magnetic 
propagation wavevector, Q=0, but with no net magne- 
tization. In neutron diffraction, the magnetic intensity 
superimposes on the nuclear Bragg peak, meaning that 
these experiments are very delicate as the magnetic in- 
tensity cx Af 2 (M is the magnetic moment) is expected to 
be very small as compared to the nuclear Bragg peaks. In 
order to detect this hidden magnetic response, polarized 
neutron experiments are then required. 



As proposed by CM. Varma [lllll^ . there are two pos- 
sibles CC phases preserving TSL. The first (the phase 
0/) has not been detected by polarized elastic neutron 
scattering experiments 0, Although it is contro- 

versial, a recent ARPES mesurement observed a dichroic 
signal in the Bi2Sr2CaCu20s+5 system consistent with 
the phase 0j/[2lj. Here, we have performed polarized 
elastic neutron scattering experiments to test the mag- 
netic moments of this second CC state, phase O//, which 
actually had never been attempted before. We success- 
fully report the first signature of a novel magnetic order 
in the pseudogap state of YBa2Cu306+a; (YBCO). The 
pattern of the observed magnetic scattering corresponds 
to the one expected in the circulating current theory of 
the pseudogap state with two current loops per Cu02 
unit-cell, phase 9// 0,^3- Alternatively, a decoration 
of the unit cell with staggered moments on the oxygen 
sites could also account for the measurements. 

All the polarized neutron diffraction measurements 
were collected on the 4FI triple-axis spectrometer at the 
Laboratoire Leon Brillouin (LLB), Saclay (France). Our 
polarized neutron diffraction setup is similar to that orig- 
inally described in [2^ with a polarized incident neutron 
with at Ei = 14.7 meV obtained with a polarizing su- 
permirror (bender) and with an Heusler analyzer (see 
also ref. 0, H^] in the context of high-Tc cuprates). 
The direction of the neutron spin polarization, P, at 
the sample position is selected by a small guide field H 
of the order of 10 G|24j]. Using that configuration, we 
monitor for each measured point the neutron scattering 
intensity in the spin-flip (SF) channel, where the mag- 
netic intensity oc M 2 is expected, and in the non-spin- 
flip (NSF) channel which measures the nuclear scatter- 
ing. To have similar counting statistics on both SF and 
NSF, we count the SF channel systematically 20 times 
longer than the NSF. We define the normalized spin-flip 
intensity as I nor m — Isf/Insf (inverse of the flipping 
ratio (FR)). With that setup, a typical flipping ratio, 
ranging between 40 and 60, is obtained. However, even 
with that high FR, the SF intensity is massively coming 
from the NSF nuclear Bragg peak through unavoidable 
polarization leakage (corresponding to about ~ 90-95% 
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TABLE I: List of samples utilized in the polarized elastic 
neutron experiments. The experiments were performed in 
the (Y,Ca)Ba2Cu3C>6+a; family for 5 samples from the un- 
derdoped (ud) to overdoped (od) part of the cuprates phase 
diagram, (t) and (d) stands for twinned and detwinned sam- 
ples, respectively. References are given where the samples 
have been described in previous neutron scattering studies. 
In contrast to the other samples, an in-plane magnetic order- 
ing occurred at Q=(l/2,l/2) in sample A with M ~ 0.05/1 b 
at 60 K[H|. 



of the SF intensity). As a very stable and homogeneous 
neutron polarization is essential through the data acqui- 
sition, all the data have been obtained in a continuous 
run versus temperature. We prove that method to be ef- 
ficient enough to see weak magnetic moments (~ 0.05/is) 
on top of nuclear Bragg peaks, see e.g. the first determi- 
nation of the A-type antiferromagnetism in Na cobaltate 
systems |25|. 

We quote the scattering wave vector as Q=(H,K,L) in 
units of the reciprocal lattice vectors, a* ~ b* = 1.63 
A -1 and c* = 0.53 A -1 . Most of the data have been 
obtained in a scattering plane where all Bragg peaks 
like Q=(0,K,L) were accessible (in twinned samples, this 
is indistinguishable from Bragg peaks with Q=(H,0,L)). 
In order to evidence small magnetic moments, measure- 
ments have been performed on the weakest nuclear Bragg 
peaks having the proper symmetry for the CC phase|26| 
(the Bragg peak Q=(0,1,1) offers the best compromise). 

We have studied 5 different samples (see Table QJ: 4 
samples in the underdoped regime and one in the over- 
doped regime. In Fig. ^a, we report the raw neutron 
intensity measured at Q=(0,1,1) for the spin flip (SF) 
channel and for the non-spin-flip (NSF) channel for an 
underdoped sample YBa 2 Cu 3 06.6(d) (sample C). The 
measurement has been done with a neutron polarization 
P//Q (sec Fig. [Jb) where the magnetic scattering is en- 
tirely spin-flip |l9ll22l E31 l24| . Between room temperature 
and a temperature T mag ~220K, the NSF and SF intensi- 
ties display the same evolution within error bars. Then, 
for T<T mag , the NSF is essentially flat whereas the SF 
intensity increases noticeably at low temperature. This 
behaviour signals the presence of a spontaneous magnetic 
order below T mQg on top of the nuclear Bragg peaks. In 
Fig. He, we show the normalized magnetic intensity as a 
function of the temperature for the 4 underdoped samples 
and the overdoped sample. For the 4 underdoped sam- 
ples, the magnetic intensity increases at low temperature 




c 

TO 





e) P-LQ 

^%^y^si. B : 



50 100 150 200 250 300 50 100 150 200 250 300 

Temperature (K) 



FIG. 1: (color online) (a): Temperature dependencies of the 
raw SF and NSF neutron intensity measured at Q=(0,1,1) 
in sample C. (b) Sketch of the scattering plane showing 
the three polarization directions discussed here, P//z cor- 
responds to the direction perpendicular to the scattering 
plane (here a*), c) Temperature dependencies of the nor- 
malized magnetic intensity, I mag , measured at Q=(0,1,1) for 
P/ /Q for the 4 underdoped samples (A,B,C,D) and the over- 
doped sample E (full points). Imag is defined as I mag {T) = 
cJnsf (300 K)[^(T) -t^(T~ S00K)] where (i) I mag 
is arbitrarily set to zero in the high temperature range (~ 
room temperature), (ii) a = 7 / I^ ) 1 aa (300K) calibrates the 
magnetic cross-sections in mbarns using the nuclear Bragg 
cross-section at Q = (0, 0, 4), 7oo| c = 7 barns. The normalized 
magnetic intensity for the Bragg peak, Q = (0,0,2), is also 
shown for samples A and C (open points), d) Temperature 
dependencies of the normalized magnetic intensity measured 
at Q=(0,1,1) (full points) (as well as Q=(0,0,2), open points) 
for P//z. e) Temperature dependencies of the normalized 
magnetic intensity, Imag, measured at Q=(0,1,1) for sample 
B for P 1 Q. 



below a certain temperature T mag whereas no magnetic 
signal is observed in the Ca-YBCO overdoped sample 
(sample E). 

We perform further measurements where the neutron 
polarization is along the complementary directions, as 
shown in Fig. Q]b, either the vertical direction P//z, or 
P _L Q but still within the horizontal scattering plane. 
The observance of the polarization selection rule for a 
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FIG. 2: a) L-scan magnetic intensity across Q=(0,1,L) in sam- 
ple A: it has been obtained using the following relation of 
measured quantities [I{SF,75K) - I(SF, 275K) j^sf^k) I 
calibrated by a (see caption Fig. 0. b) Temperature de- 



pendencies of the magnetic intensity, 
peaks L=0,l,2 in sample A. 



for various Bragg 



magnetic signal, Ip//Q = Ip//z + ^P-LQj m the three po- 
larizations, as shown in Fig. ^c^d and^ e f° r sample 
B, unambiguously demonstrates the magnetic origin of 
the low temperature signal. More precisely, in the P//z, 
configuration, only magnetic moments within the hori- 
zontal scattering plane but sti ll perp endicular to Q are 
observed in SF channel |H HI Hi Hj . For Q = (0, 1, 1), 
this means that we mostly probe the magnetic moments 
parallel to the c* axis. In the 4 underdoped samples, we 
observe a similar onset of the magnetic order below T mag 
for P//z (Fig. G]d). This demonstrates that the deduced 
magnetic moment has a well-defined component perpen- 
dicular to the CuC>2 plane. However, a closer comparison 
with both polarizations reveals that their intensities do 
not simply match. This underlines that the magnetic 
moment also exhibits an in-plane component (within the 
CuC>2 plane) as the cross-section in Fig. ^ c is larger than 
the one in Fig. ^d. Combining all measured polariza- 
tions in the different samples, one can then estimate a 
mean angle between the direction of the moments with 
the c* axis to be 4> — 45° ± 20° valid for all samples. 

As shown in Fig. \I\c, the typical cross-section of the 
magnetic order is ~ 1-2 mbarns, i.e. ~ 1CP 4 of the 
strongest Bragg peaks. This explains why such a mag- 
netic order was not reported before with unpolarized neu- 
tron diffraction. Due to these experimental limitations, 
we do not perform a detailed and quantitative determi- 
nation of magnetic structure for which further work is 
needed. However, some qualitative aspects can be briefly 
discussed. First, we perform a scan along the L-direction 
in the SF channel across the Bragg peak (Fig. [21a) where 
the difference in temperature between T=75 K and 275 
K has been taken to remove the effect of the polariza- 
tion leakage. The observed magnetic peak is resolution 
limited, showing that the magnetic order is character- 



FIG. 3: (color online) a) Cuprate superconductors phase di- 
agram as a function of hole doping, n/,, deduced from the 
SC temperature using the empirical relation Tc/T™ ax = 
1 — 82.6(rih — 0.16) 2 [3l|. The white points show Tmag (see 
table rU- The color map shows the quantity 8R(T) = 1 — 
[p a b{T) — Pab(0)]/[cxT) deduced from the resistivity measure- 
ments in YBCO|g: the change of colors indicates the depar- 
ture from the T-linear behaviour, SR(T) — represented in 
blue, at high temperature. SR(T) 7^ defines the pseudo- 
gap state, b) Circulating current phase, On, in the CUO2 
in t he p s 

superconducting cuprates[Tl1.ll2|. 
TSL. 



A spin model preserving 



ized by long range 3D correlations at T=75 K. Second, 
by looking at other Bragg peaks along c* (Fig. |21b), we 
found that the magnetic intensity is not uniformly dis- 
tributed versus L, meaning that i) the magnetic intensity 
does not arise from the Cu-0 chains, and ii) the moments 
arrangement within a bilayer appears to be mainly paral- 
lel. This directly arises from the hierachy of the observed 
magnetic intensities (intensity at L=0 is larger than at 
L=2, Fig. 0b). Finally, using the observed magnetic 
cross-section (Fig. ^c) and a weakly momentum depen- 
dent form factor, one can deduce a typical magnitude of 
ordered magnetic moment of M ~ 0.05 to 0.1//B with 
the moment decreasing with increasing doping in the 4 
samples. 

Therefore, we observe an unusual magnetic order in 
a temperature and doping range that cover the range 
where the pseudogap state is observed in YBCO. Our 
data do not contradict previous unsuccessfull polarized 
neutron reports^^ as the Bragg spots where the effect 
is observed are along a direction at 45° from the one 
previously studied|2y|. The deduced T mag , defined as 
the change of slope in the normalized intensity I mag , de- 
creases with increasing doping (see table Pi . It matches 
the pseudogap temperature, T*, of the resistivity data 
in YBCO |(j as shown on Fig. [21 The occurrence of a 
magnetic order in this temperature and doping ranges 
points towards a magnetic signature of an hidden order 
parameter associated with the pseudogap state. As all 
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anomalous physical properties evidencing the pseudogap, 
the temperature dependence of the magnetic order does 
not exhibit a marked change at T mag . Being on top of 
nuclear Bragg peaks, that magnetic order does not break 
TSL, indicating a zero magnetic propagation wavevec- 
tor, Q = 0. As shown in Fig. QJ ; no magnetic intensity 
occurs below T mag at the Bragg peak Q = (0,0,2), rul- 
ing out a ferromagnetic order. The absence of breaking 
of TSL points towards a magnetic pattern of antipar- 
allel magnetic moments within each unit cell. Among 
the proposed order parameters, only one gives magnetic 
scattering at Q = (10L)[= (OIL)]: it is the orbital mo- 
ments arising from the circulating current phase with 2 
current loops per CuC>2 unit cell, 6// (Fig. 0b) flilllij]. 
Another possibility could be a model with colinear spin 
moments located at the oxygen site ferromagnetic along 
both directions of the square lattice but antiparallcl each 
other as sketched on Fig. |3|c. Any other model charac- 
terized by a decoration of the CuC>2 plaquette would also 
give rise to a magnetic contribution at the proper Bragg 
spots. 

From our present measurements, one cannot distin- 
guish between these two models. Only a detailed study 
of magnetic form factors would allow to differentiate the 
scattering from spin and orbital moments. However, 
some arguments can be given from the observed moments 
direction, which is not within the CuC>2 plane neither 
perpendicular to it. Clearly, an in-plane magnetic com- 
ponent is not expected within the orbital moments pic- 
ture of currents flowing within perfectly flat CuC>2 planes. 
However, due to the dimpling of Cu02 planes in YBCO, 
the moments can be tilted by about 11° from the c axis 
as the effective moments at the centers of the O-Cu-0 
plaquettes are perpendicular to these plaquettes. Within 
an orbital moment picture, spin degree of freedom might 
also play a role in producing in-plane magnetic moments, 
for instance, by spin-orbit scattering or in relation to 
chiral spin states associated with flux phases [28|. Al- 
ternatively, if considering spin models, one would rather 
expect moments lying within the Cu02 plaquette as it 
is the case for copper spins in undoped cuprates. A rea- 
son should be found to explain why the moments ex- 
hibit an out-of-plane component. Whatever the origin 
of the observed order, its pattern challenges the single 
band Hubbard picture commonly used to describe high- 
Tc cuprates. At very least, oxygen orbitals need to be 
included to determine the minimal effective Hamiltonian 
for the cuprates. 

In conclusion, we report a first signature of an unusual 
magnetic order in several YBa2Cu306+j; samples match- 
ing the pseudogap behaviour in underdoped cuprates 
(Fig. |2J). Such an observation points towards the ex- 
istence of an hidden order parameter for the pseudogap 
phase in high-Tc superconductors. Importantly, our ex- 
periment reveals that a 3D long range order does not 
break the translational symmetry of the lattice and im- 



plies a decoration of the unit cell with staggered spin or 
orbital moments. The symmetry of the observed order 
corresponds to the one expected in orbital moments em- 
anating from a circulating current state (ill fl2j | . 
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